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ABSTRACT 

We analyse the large-scale correlation function of the 6dF Galaxy Survey (6dFGS) 
and detect a Baryon Acoustic Oscillation (BAO) signal. The 6dFGS BAO detection 
allows us to constrain the distance-redshift relation at z e s — 0.106. Wc achieve a 
distance measure of Dv(z c s) = 456 ±27 Mpc and a measurement of the distance ratio, 
r s (zd)/Dv{z c ff) = 0.336 ± 0.015 (4.5% precision), where r s {zd) is the sound horizon 
at the drag epoch Zd- The low effective redshift of 6dFGS makes it a competitive and 
independent alternative to Cepheids and low- 2; supernovac in constraining the Hubble 
constant. We find a Hubble constant of Ho = 67 ±3.2 km s _1 Mpc -1 (4.8% precision) 
that depends only on the WMAP-7 calibration of the sound horizon and on the galaxy 
clustering in 6dFGS. Compared to earlier BAO studies at higher redshift, our analysis 
is less dependent on other cosmological parameters. The sensitivity to Hq can be used 
to break the degeneracy between the dark energy equation of state parameter w and 
H in the CMB data. We determi ne that w = -0 97 ± 0.13, using only WMAP-7 and 
BAO data from both 6dFGS and lPercival et all (l2010f ). 

We also discuss predictions for the large scale correlation function of two fu- 
ture wide-angle surveys: the WALLABY blind HI survey (with the Australian SKA 
Pathfinder, ASKAP), and the proposed TAIPAN all-southern-sky optical galaxy sur- 
vey with the UK Schmidt Telescope (UKST). We find that both surveys are very likely 
to yield detections of the BAO peak, making WALLABY the first radio galaxy survey 
to do so. We also predict that TAIPAN has the potential to constrain the Hubble 
constant with 3% precision. 

Key words: surveys, cosmology: observations, dark energy, distance scale, large scale 
structure of Universe 



1 INTRODUCTION 

The current standard cosmological model, ACDM, assumes 
that the initial fluctuations in the distribution of matter 
were seeded by quantum fluctuations pushed to cosmologi- 
cal scales by inflation. Directly after inflation, the universe 
is radiation dominated and the baryonic matter is ionised 
and coupled to radiation through Thomson scattering. The 
radiation pressure drives sound-waves o riginating from over- 
densities in the matter distribution (|Peebles fc Yul Il970l : 
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ISunvaev & Zeldovicbl |l97Ct iBond fc Efstathioul 1 19871 ). At 
the time of recombination (2* « 1090) the photons decouple 
from the baryons and shortly after that (at the baryon drag 
epoch Zd ~ 1020) the sound wave stalls. Through this pro- 
cess each over-density of the original density perturbation 
field has evolved to become a cen trally peaked perturbation 
surrounded by a spherical shell (Bashinsky fc Bertschingeil 
120011 120021 ; lEisenstein. Seo fc White! 120071 ). The radius of 
these shells is called the sound horizon r 3 . Both over-dense 
regions attract baryons and dark matter and will be pre- 
ferred regions of galaxy formation. This process can equiv- 
alently be described in Fourier space, where during the 
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photon-baryon coupling phase, the amplitude of the baryon 
perturbations cannot grow and instead undergo harmonic 
motion leading to an oscillation pattern in the power spec- 
trum. 

After the time of recombination, the mean free path 
of photons increases and becomes larger than the Hubble 
distance. Hence from now on the radiation remains almost 
undisturbed, eventually becoming the Cosmic Microwave 
Background (CMB). 

The CMB is a powerful probe of cosmology due to the 
good theoretical understanding of the physical processes de- 
scribed above. The size of the sound horizon depends (to 
first order) only on the sound speed in the early universe 
and the age of the Universe at recombination, both set 
by th e physical matter and baryon densities, Q. m h 2 and 
Q b h 2 (|Eisenstein fc Hull 19981 1. Hence, measuring the sound 
horizon in the CM B gives extremely ac curate constraints 
on these quantities jKomatsu et aill2010ft . Measurements of 
other cosmological pa rameters often show deg eneracies in 
the CMB data alone (|Efstathiou fc Bondlll999ft . especially 
in models with extra parameters beyond flat ACDM. Com- 
bining low redshift data with the CMB can break these de- 
generacies. 

Within galaxy redshift surveys we can use the correla- 
tion function, £, to quantify the clustering on different scales. 
The sound horizon marks a preferred separation of galaxies 
and hence predicts a peak in the correlation function at the 
corresponding scale. The expected enhancement at s = r 3 
is only A£ « 1CT 3 6 2 (1 + 2/3/3 + (3 2 /5) in the galaxy cor- 
relation function, where b is the galaxy bias compared to 
the matter correlation function and /3 accounts for linear 
redshift space distortions. Since the signal appears at very 
large scales, it is necessary to probe a large volume of the 
universe to decrease s ample variance, which dominates the 
error on these scales ( Tegmarkl]l997| ; iGoldberg fc Strauss! 
1 19981 : lEisenstein. Hu fe Tegmarklll99Sft . 

Very interesting for cosmology is the idea 
of usi ng the sound horizon sc ale a s a standard 
ruler (lEisenstein. Hu fc Tegmark 19981; Coorav et al.l 



same data found similar r e sults using the correlation func- 
tion Z(s) jMartinezI 120091; iGaztanaga Cabre fc Huil 120091 ; 



120011 : ISeo fc Eisensteinll2003l : lBlake fc Glazebrookll2003ft . A 
standard ruler is a feature whose absolute size is known. By 
measuring its apparent size, one can determine its distance 
from the observer. The BAO signal can be measured in 
the matter distribution at low redshift, with the CMB cal- 
ibrating the absolute size, and h ence the distance-redshift 
relation can be mapped (see e.g. iBassett fc Hloze"k1 (|2009ft 

for a summary). 

The Sloan Digital Sky Survey (SDSS; lYork et all 
2000ft, and the 2d F Galaxy Redshift Survey (2dFGRS; 



Colless et al.1 l200lft were the first redshift surveys which 



have directly detected the BAO signal. Recently the Wig- 
gleZ Dark Energy S urvey has reporte d a BAO measurement 



glez; JJark Jincrgy Survey lias reported 
at r edshift z = 0.6 (iBlake et ablhoilft . 

lEisenstein et al. I l|2005ft were able to constrain the 
distance-redshift relation to 5% accuracy at an effective 
redshift of z c g — 0.35 using an early data release of the 
SDSS-LRG sample containing « 47 000 galaxies. Subse- 
quent studies using the final SDSS-LRG sample and com- 
bining it with the SDSS-main and the 2dFGRS sample 
were able to improve on this measurement and constrain 
the distance-redsh ift relation at z e g = 0.2 and 2 e fi = 0.35 
with 3% accuracy (|Percival et aillioioft . Other studies of the 



lLabini et a.1.1 12009| ; ISanchez et al.l 120091; iKazin et al.l 20 10ft 

the power spectrum P(k) dCole et al. 120051 ; Tegmark et al.l 
120061 ; iHuetsil 120061 ; iReid et ail l201Cft . the projected cor- 
relation function w(r v ) of photometric red shift samples 
jPadmanabhan et al.ll2007l ; IBlake et aill2007ft and a cluster 
sample based on the SDSS photo metric data dHuetsi 20091) . 
Severa l years earlier a study by iMiller. Nichol fc Batuskl 
|200lft found first hints of the BAO feature in a combination 
of smaller datasets. 

Low redshift distance measurements can directly mea- 
sure the Hubble constant Ho with a relatively weak depen- 
dence on other cosmological parameters such as the dark en- 
ergy equation of state parameter w. The 6dF Galaxy Survey 
is the biggest galaxy survey in the local universe, covering 
almost half the sky. If 6dFGS could be used to constrain 
the redshift-distance relation through baryon acoustic oscil- 
lations, such a measurement could directly determine the 
Hubble constant, depending only on the calibration of the 
sound horizon through the matter and baryon density. The 
objective of the present paper is to measure the two-point 
correlation function on large scales for the 6dF Galaxy Sur- 
vey and extract the BAO signal. 

Many cosmological parameter studies add a prior on 
Ho to help break degeneracies. The 6dFGS derivation of Ho 
can provide an alternative source of that prior. The 6dFGS 
Ho-measurement can also be used as a consistency check of 
other low redshift distance calibrators such as Cepheid vari- 
ables and Type la supern ovae (through the so ca l led distance 
ladde r technique; see e.g. iFreedman et al.ll200ll ; iRiess et ail 
l201lft . Compared to these more classical probes of the Hub- 
ble constant, the BAO analysis has an advantage of simplic- 
ity, depending only on Q. m h 2 and Qbh 2 from the CMB and 
the sound horizon measurement in the correlation function, 
with small systematic uncertainties. 

Another motivation for our study is that the SDSS data 
after data release 3 (DR3) show more correlation on large 
scales than expected by ACDM and have no sign of a cross- 
over to negative f up to 200/t~ 1 Mpc (the ACDM prediction 
is 140/2T 1 Mpc) |Kazin et ail 12010ft . It could be that the 
LRG sample is a rather unusual realisation, and the addi- 
tional power just reflects sample variance. It is interesting 
to test the location of the cross-over scale in another red 
galaxy sample at a different redshift. 

This paper is organised as follows. In Section [2] we 
introduce the 6dFGS survey and the A'-band selected 
sub-sample used in this analysis. In Section [3] we explain 
the technique we apply to derive the correlation function 
and summarise our error estimate, which is based on 
log-normal realisations. In Section 2] we discuss the need 
for wide angle corrections and several linear and non-linear 
effects which influence our measurement. Based on this 
discussion we introduce our correlation function model. In 
Section [5] we fit the data and derive the distance estimate 
Dv(z e s). In Section [6] we derive the Hubble constant and 
constraints on dark energy. In Section [7] we discuss the 
significance of the BAO detection of 6dFGS. In Section [8] 
we give a short overview of future all-sky surveys and their 
power to measure the Hubble constant. We conclude and 
summarise our results in Section [9] 
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Throughout the paper, we use r to denote real space 
separations and s to denote separations in redshift space. 
Our fiducial model assume a flat universe with Q,™ = 0.27, 
w Rd = -1 and fi£ d = 0. The Hubble constant is set to 
Ho = 100ft km s -1 Mpc -1 , with our fiducial model using 
ft fld = 0.7. 



2 THE 6DF GALAXY SURVEY 

2.1 Targets and Selection Function 

The galaxies used in this analysis were selected to K ^ 
12.9 from the 2MASS E xtended Source Catalog (2MASS 
XSC; IJarrett et afll200Ch and combined with redshift data 
from the 6dF Galaxy Survey (6dFGS; I Jones et al.l l2009h . 
The 6dF Galaxy Survey is a combined redshift and peculiar 
velocity survey covering nearly the entire southern sky with 
\b\ < 10°. It was undertaken with the Six-Degree Field (6dF) 
multi-fibre instrument on the UK Schmidt Telescope from 
2001 to 2006. The median redshift of the survey is z = 0.052 
and the 25% : 50% : 75% perce ntile comple t eness values 
are .61 : 0.79 : 0.92. Papers by Ijones etahl (|2004 l20od . 
120091 ) describe 6dFGS in full detail, including comparisons 
between 6dFGS, 2dFGRS and SDSS. 

Galaxies were excluded from our sample if they resided 
in sky regions with completeness lower than 60 percent. Af- 
ter applying these cuts our sample contains 75 117 galaxies. 
The selection function was derived by scaling the survey 
completeness as a function of magnitude to match the inte- 
grated on-sky completeness, us ing mean galaxy cou nts. This 
method is the same adopt ed bvlColless et al] |200ll ) for 2dF- 
GRS and is explained in Ijones et al.l ( 20061 ) in detail. The 
redshift of each object was checked visually and care was 
taken to exclude foreground Galactic sources. The derived 
completeness function was used in the real galaxy catalogue 
to weight each galaxy by its inverse completeness. The com- 
pleteness function was also applied to the mock galaxy cat- 
alogues to mimic the selection characteristics of the survey. 
Jones et al. (in preparation) describe the derivation of the 
6dFGS selection function, and interested readers are referred 
to this paper for a more comprehensive treatment. 



Az = 0.001 




Figure 1. Redshift distribution of the data (black solid line) and 
the random catalogue (black dashed line) . The weighted distribu- 
tion (using weights from eq. [6]l is shifted to higher redshift and has 
increased shot noise but a smaller error due to sample variance 
(blue solid and dashed lines). 



The volume of the 6dF Galaxy Survey is approxi- 
mately as large as the volume covered by the 2dF Galaxy 
Reds hift Survey, with a sample density similar to SDSS- 
DR7 (jAbazaiian et al]|2009h . |Percival et all (|2010l ) reported 
successful BAO detections in several samples obtained from 
a combination of SDSS DR7, SDSS-LRG and 2dFGRS 
with effective volumes in the range 0.15 — 0.45ft -3 Gpc 3 
(using P = 10 000ft -3 Mpc 3 ), while the original detec- 
tion by lEisenstein et al.l (|2005l ) used a sample with V c s = 



0.38ft -3 Gpc 3 (using P = 40 000ft 



3 Mpc 3 ^ 



3 CLUSTERING MEASUREMENT 

We focus our analysis on the two-point correlation func- 
tion. In the following sub-sections we introduce the tech- 
nique used to estimate the correlation function and outline 
the method of log-normal realisations, which we employed 
to derive a covariance matrix for our measurement. 



2.2 Survey volume 

We calcula ted the ef f ective volume of the survey using the 
estimate of lTegmarkl (|l997T ) 



VeB = 



n(x)P 



1 + n(x)P Q 



(1) 



where n{x) is the mean galaxy density at position x, de- 
termined from the data, and Po is the characteristic power 
spectrum amplitude of the BAO signal. The parameter Po 
is crucial for the weighting scheme introduced later. We find 
that the value of Po = 40 000ft -3 Mpc 3 (corresponding to 
the value of the galaxy power spectrum at k « 0.06ft Mpc -1 
in 6dFGS) minimises the error of the correlation function 
near the BAO peak. 

Using Po = 40 000ft -3 Mpc 3 yields an effective volume 
while using instead Po = 10 000ft -3 Mpc 3 



of 0.08ft -3 Gpc 3 , 



(corresponding to k ~ 0.15ft Mpc ) gives an effective vol- 
ume of 0.045ft -3 Gpc 3 . 



3.1 Random catalogues 

To calculate the correlation function we need a random sam- 
ple of galaxies which follows the same angular and redshift 
selection function as the 6dFGS sample. We base our ran- 
dom catalogue generation on the 6dFGS luminosity func- 
tion of Jones et al. (in preparation), where we use ran- 
dom numbers to pick volume-weighted redshifts and lumi- 
nosity function-weighted absolute magnitudes. We then test 
whether the redshift-magnitude combination falls within the 
6dFGS A'-band faint and bright apparent magnitude limits 
(8.75 < K < 12.9). 

Figure [T] shows the redshift distribution of the 6dFGS 
A'-selected sample (black solid line) compared to a random 
catalogue with the same number of galaxies (black dashed 
line). The random catalogue is a good description of the 
6dFGS redshift distribution in both the weighted and un- 
weighted case. 
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where the ratio n r /rid is given by 



0.05^ 
0.04 r 
0.03 r 
0.02^ 
0.01 

0^ 
-0.01 - 



6dFGS data 
best fit 
I2 m h 2 = 0.12 
I2 m h 2 = 0.15 
no-baryon fit 




-0.02 



20 



60 



80 100 120 
s [h 1 Mpc] 



140 160 180 200 



Figure 2. The large scale correlation function of 6dFGS. The 
best fit model is shown by the black line with the best fit value of 
tt m h 2 = 0.138 ± 0.020. Models with different n m h 2 are shown by 
the green line {Q.mh 2 = 0.12) and the blue line (Q. m h 2 = 0.15). 
The red dashed line is a linear CDM model with Qi,h 2 = (and 
Ci m h 2 = 0.1), while all o ther models use the WMAP-7 best fit 
value of n b h 2 = 0.02227 jKomatsu et alj|2010h . The significance 
of the BAO detection in the black line relative to the red dashed 
line is 1Aa (see Section 0. The error-bars at the data points 
are the diagonal elements of the covariance matrix derived using 
log-normal mock catalogues. 



3.2 The correlation function 

We turn the measured redshift into co-moving distance via 

dz 



Dc(z) = -- 

H o Jo 



E(z') 



with 



E(z) = [Sl™(l + z) 3 + fi£ d (l + z ) 2 



+ nr(i + z) 



3(lW ld )y|l/2 



(2) 



(3) 



where the curvature fl^. d is set to zero, the dark energy den- 



sity is given by f^ d 



1 — njjf and the equation of state for 
dark energy is w na = — 1. Because of the very low redshift 
of 6dFGS, our data are not very sensitive to fife, w or any 
other higher dimensional parameter which influences the ex- 
pansion history of the universe. We will discuss this further 
in Section \5. 31 

Now we measure the separation between all galaxy pairs 
in our survey and count the number of such pairs in each sep- 
aration bin. We do this for the 6dFGS data catalogue, a ran- 
dom catalogue with the same selection function and a com- 
bination of data-random pairs. We call the pair-separation 
distributions obtained from this analysis DD(s), RR(s) and 
DR(s), respectively. The binning is chosen to be from 
10/iT 1 Mpc up to 190/T 1 Mpc, in lO/i" 1 Mpc steps. In the 
analysis we used 30 random catalogues with the same size as 
the data c atalogue. The redshift correlation function itself 
is given bv lLandv fc Szalavl (| 19931 ): 



Cdata(s) = 1 + 



DD(s) 
RR(s) \nd 



RR(s) 



^ , (4) 
rid ! 



rid 



(5) 



and the sums go over all random (N r ) and data 
( Nd) galaxies. We use the inver se density weighting 
of iFeldman, Kaiser fc Peacock! (119941 ): 

m{g) = i+Sw (6) 

with P = 40 000/i 3 Mpc" 3 and d being the inverse com- 
pleteness weighting for 6dFGS (see Section ^. H and Jones et 
al., in preparation). This weighting is designed to minimise 
the error on the BAO measurement, and since our sample 
is strongly limited by sample variance on large scales this 
weighting results in a significant improvement to the analy- 
sis. The effect of the weighting on the redshift distribution 
is illustrated in Figure [T] 

Other authors have used the so called ^-weighting 
which optimises the error over all scales by weighting eac h 
scale differently (e.g. lEfstathioulll988l : lLovedav et al.lll995l ). 
In a magnitude limited sample there is a correlation be- 
tween luminosity and redshift, which establishes a correla- 
tion between bias and redshift (|Zehavi et al.ll2005l ). A scale- 
dependent weighting would imply a different effective red- 
shift for each scale, causing a scale dependent bias. 

Finally we considered a luminosity dep e ndent weighting 
as suggested bv lPercival. Verde fc Peac ock (2004|). However 
the same authors found that explicitly accounting for the 
luminosity-redshift relation has a negligible effect for 2dF- 
GRS. We found that the effect to the 6dFGS correlation 
function is la for all bins. Hence the static weighting of 
eq. [6] is sufficient for our dataset. 

We also include an integral constraint correction in the 
form of 



Cdata(s) = Cdata(s) + »C, 



where ic is defined as 



J2 S RR(s)£model(s) 

E s RR{s) 



(7) 



(8) 



The function RR(s) is calculated from our mock catalogue 
and £ mo dei(s) is a correlation function model. Since ic de- 
pends on the model of the correlation function we have to 
re-calculate it at each step during the fitting procedure. 
However we note that ic has no significant impact to the 
final result. 

Figure [2] shows the correlation function of 6dFGS at 
large scales. The BAO peak at « 105/i _1 Mpc is clearly 
visible. The plot includes model predictions of different cos- 
mological parameter sets. We will discuss these models in 
Section El 



3.3 Log-normal error estimate 

To obtain reliable error-bar s for the correlation f unction we 
use log-normal realisation s (| Coles fc Jonedli~99ll : I Cole et al.l 
120051 : iKitaura et al.l |2009| ) . In what follows we summarise 
the main steps, but refer the interested reader to Ap- 
pendix [X] in which we give a detailed explanation of how we 
generate the log-normal mock catalogues. In Appendix [B] 
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Figure 3. Correlation matrix derived from a covariance matrix 
calculated from 200 log-normal realisations. 



we compare the log-normal errors with jack-knife estimates. 



N — 1 



(11) 



Here, £n(*i) is the correlation function estimate at separa- 
tion Si and the sum goes over all iV log-normal realisations. 
The mean value is defined as 



l N 



(12) 



The case i — j gives the error (ignoring correlations between 
bins, of = Ca). In the following we will use this uncertainty 
in all diagrams, while the fitting procedures use the full co- 
variance matrix. 

The distribution of recovered correlation functions in- 
cludes the effects of sample variance and shot noise. Non- 
linearities are also approximately included since the distri- 
bution of over-densities is skewed. 

In Figure |3] we show the log-normal correlation matrix 
nj calculated from the covariance matrix. The correlation 
matrix is defined as 



\J CiiCjj 



(13) 



where C is the covariance matrix (for a comparison to jack- 
knife errors see appendix |B|) . 



Log-normal realisations of a galaxy survey are usually 
obtained by deriving a density field from a model power 
spectrum, P(k), assuming Gaussian fluctuations. This den- 
sity field is then Poisson sampled, taking into account the 
window function and the total number of galaxies. The as- 
sumption that the input power spectrum has Gaussian fluc- 
tuations can only be used in a model for a density field 
with over-densities <C 1. As soon as we start to deal with 
finite rms fluctuations, the Gaussian model assigns a non- 
zero probability to regions of negative density. A log-normal 
random field LN(x), can avoid this unphysical behaviour. 
It is obtained from a Gaussian field G(x) by 



LN(x) = exp[G(f) 



(9) 



which is positive-definite but approaches 1 + G(x) when- 
ever the perturbations are small (e.g. at large scales). Cal- 
culating the power spectrum of a Poisson sampled density 
field with such a distribution will reproduce the input power 
spectrum convolved with the window function. As an input 
power spectrum for the log-normal field we use 



Pm(fc) = AP lin (fc)exp[-(fc/fc,r 



(10) 



where A = b 2 (l + 2/3/3 + /3 2 /5) accounts for the linear 
bias and the linear redshift space distortions. Pu n (k) is a 
linear m odel power spectru m in real space obtained from 
CAMB (|Lewis et all boOut ) and P n i(k) is the non-linear 
power spectrum in redshift space. Comparing the model 
above with the 6dFGS data gives A — 4. The damping pa- 
rameter fc* is set to fc* = 0.33/t Mpc -1 , as found in 6dFGS 
(see fitting results later). How well this input model matches 
the 6dFGS data can be seen in Figure [9] 

We produce 200 such realisations and calculate the cor- 
relation function for each of them, deriving a covariance ma- 



4 MODELLING THE BAO SIGNAL 

In this section we will discuss wide-angle effects and non- 
linearities. We also introduce a model for the large scale 
correlation function, which we later use to fit our data. 



4.1 Wide angle formalism 

Th e mode l of lin ear redshift space distortions introduced 
by iKaiserl (|l987h is based on the plane parallel approxi- 
mation. Earlier surveys such as SDSS and 2dFGRS are at 
sufficiently high redshift that the maximum opening an- 
gle between a galaxy pair remains small enough to en- 
sure the plane parallel approximation is valid. However, the 
6dF Galaxy Survey has a maximum opening angle of 180° 
and a lower mean redshift of z fa 0.1 (for our weighted 
sample) and so it is necessary to test the validity of the 
plane parallel approximation. The wide angle description 
of reds hift space distortions has been laid out in severa l 
papers dSzalav et af]|l997l: ISzapudil 12004] ; lMatsubarall2004l ; 
IPapai fc Szapudil 120081 ; iRaccanelli et alj l201Ch . which we 
summarise in Appendix [Cl 

We find that the wide-angle corrections have only a very 
minor effect on our sample. For our fiducial model we found a 
correction of A£ = 4- 10~ 4 in amplitude at s = 100ft -1 Mpc 



and A£ = 4.5 ■ 10~ 4 at s = 200/i _1 Mpc, (Figure [C2] in 
the appendix). This is much smaller than the error bars on 
these scales. Despite the small size of the effect, we never- 
theless include all first order correction terms in our correla- 
tion function model. It is important to note that wide angle 
corrections affect the correlation function amplitude only 
and do not cause any shift in the BAO scale. The effect of 
the wide-angle correction on the unweighted sample is much 
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greater and is already noticeable on scales of 20/i _1 Mpc. 
Weighting to higher redshifts mitigates the effect because it 
reduces the average opening angle between galaxy pairs, by 
giving less weight to wide angle pairs (on average). 

4.2 Non-linear effects 

There are a number of non-linear effects which can poten- 
tially influence a measurement of the BAO signal. These in- 
clude scale-dependent bias, the non-linear growth of struc- 
ture on smaller scales, and redshift space distortions. We 
discuss each of these in the context of our 6dFGS sample. 

As the universe evolves, the acoustic signature in the 
correlation function is broadened by non-linear gravitational 
structure formation. Equivalently we can say that the higher 
harmonics in the power spectrum, which represe nt smaller 
scales, are erased (|Eisenstein. Seo fc White]|2007r ). 

The early universe physics, which we discussed briefly in 
the introduction, is well understood and several authors have 
produced software packages (e.g. C MBFAST and CAMB) 
and published fitting functions (e.g lEisenstein fe Hrj|l998h 
to make predictions for the correlation function and power 
spectrum using thermodynamical models of the early uni- 
verse. These models already include the basic linear physics 
most relevant for the BAO peak, fn our analy sis we use 
the CAMB software package (|Lewis et al.ll200fj ). The non- 
linear evolution of the power sp ectrum in CAMB is cal- 
culated using the hale-fit code (|Smith et all I2003T I. This 
code is calibrated by n-body simulations and can describe 
non-linear effects in the shape of the matter power spec- 
trum for pure CDM models to an accuracy of around 
5 - 10% jHeitmann et al.l |2010» . However, it has previ- 
ously been shown that this non-linear model is a poor 
description of the non-linear effects around the BAO 



peak (|Crocce fc Scoccimarr ol l200gf ). We therefore decided to 
use the linear model output from CAMB and incorporate 
the non-linear effects separately. 

All non-linear effects influencing the correlation func- 
tion can be approximated by a con volution with a Gaus- 
sian damping factor expl — (rk* /2) 2 ] (|Eisenstein et al]|2007l ; 
lEisenstein. Seo fc White! 120071 ). where fc, is the damping 
scale. We will use this factor in our correlation function 
model introduced in the next section. The convolution with 
a Gaussian causes a shift of the peak position to larger scales, 
since the correlation function is not symmetric around the 
peak. However this shift is usually very small. 

All of the non-linear processes discussed so far are not 
at the fundamental scale of 105/i _1 Mpc but are instead at 
the cluster-formation scale of up to 10h~ Mpc. The scale of 
105h~ 1 Mpc is far larger than any known non-linear effect in 
cosmology. This has led some authors to the conclusion that 
the peak will not be shift ed significantly, but rath e r only 
blurred out. For example, lEisenstein. Seo fc White! (|2007T ) 
have argued that any systematic shift of the acoustic scale 
in real space must be small (< 0.5%), even at z = 0. 



to 1% 



Ho wever, several authors re p ort possible shifts of up 
(iGuzik fc Bernsteinl 120071; ISmith et all |200S| , 120071 ; 



lAngulo et al.l 120081 ). ICrocce fc Scoccimar ro (20o|)used 



normalised perturbation theory (RPT) and found percent- 
level shifts in the BAO peak. In addition to non-linear evo- 
lution, they found that mode-coupling generates additional 
oscillations in the power spectrum, which are out of phase 



with the BAO oscillations predicted by linear theory. This 
leads to shifts in the scale of oscillation nodes with re- 
spect to a smooth spectrum. In real space this corresponds 
to a peak shift towards sma l ler sc ales. Based on their re- 
sults, [Crocce^&^Scoccimarr^ (|2008l ) propose a model to be 
used for the correlation function analysis at large scales. We 
will introduce this model in the next section. 



4.3 Large-scale correlation function 

To model the correlation f unctio n o n large scales, we fol - 
low ICrocce fc Scoccimarrcl (|2008t ) and lSanchez et al.l (|2008T ) 
and adopt the following parametrisatiorfjj 



£modol(s) = B(s)b 2 



£(s)*G(r)+£l(r) 



ae( a ) 



ds 



(14) 



Here, we decouple the scale dependency of the bias B(s) 
and the linear bias b. G(r) is a Gaussian damping term, ac- 
counting for non-linear suppression of the BAO signal. £(s) is 
the linear correlation function (including wide angle descrip- 
tion of redshift space distortions; eq. IC4I in the appendix) . 
The second term in eq [2] accounts for the mode-coupling of 
different Fourier modes. It contains d£(s)/ds, which is the 
first derivative of the redshift space correlation function, and 
£i( r )i which is defined as 

£i>) = ^2_/ o dk kPimWMrk), (15) 

with ji(x) be i ng the spherical Bessel function of order 1. 
ISanchez et al.l (|2008h used an additional parameter Amc 
which multiplies the mode coupling term in equation 1141 
We found that our data is not good enough to constrain 
this paramete r, and hence adopted Amc = 1 as in the orig- 
inal model bv lCrocce fc Scoccimarrcl (|2008l ). 

In practice we generate linear model power spectra 
fiin(fc) from CAMB and convert them into a correlation 
function using a Hankel transform 



i 

2^2 



dk fc Pii n (fc)jo(rfc), 



(16) 



where jo(x) = sin(x)/x is the spherical Bessel function of 
order 0. 

The *-symbol in eq. [14] is only equivalent to a convo- 
lution in the case of a 3D correlation function, where we 
have the Fourier theorem relating the 3D power spectrum 
to the correlation function. In case of the spherically aver- 
aged quantities this is not true. Hence, the *-symbol in our 
equation stands for the multiplication of the power spectrum 
with G(k) before transforming it into a correlation function. 
G(fc) is defined as 



G(fc) = exp [-(fc/fc,) 2 ] 



(17) 



with the property 

G(fc) as fc -> oo. (18) 
The damping scale fc, can be calculated from linear the- 



1 note that r = s, the different letters just specify whether the 
function is evaluated in redshift space or real space. 
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Table 1. This table contains all parameter constraints from 
6dFGS obtained in this paper. The priors used to derive these 
parameters are listed in square brackets. All parameters as- 
sume n b h 2 = 0.02227 and in cases where a prior on ^If^h is 
used, we adopt the WMA P-7 Markov chain probability distribu- 
tion llKomats u et al.||201C ). A( z r ff) is the acoustic parameter de- 
fined byjEisensteh^^Lf 20051 ) (see equation [27] in the text) and 
R(z e g) is the distance ratio of the 6dFGS BAO measurement to 
the last-scattering surface. The most sensible value for cosmolog- 
ical parameter constraints is r s {z^)/ Dy{z e fi ), since this measure- 
ment is uncorrelated with Q, m h? . The effective redshift of 6dFGS 
is z e ff = 0.106 and the fitting range is from 10 — 190/i — 1 Mpc. 

Summary of parameter constraints from 6dFGS 



n m h 2 

Dv(zctt) 
r s (z d )/D v (z eff ) 

R(Zes) 
A(z eff ) 



0.138 ±0.020 (14.5%) 
456 ± 27 Mpc (5.9%) 
459 ± 18 Mpc (3.9%) 
0.336 ±0.015 (4.5%) 
0.0324 ± 0.0015 (4.6%) 
0.526 ±0.028 (5.3%) 



[tor, 



prior] 



fl m 0.296 ± 0.028 (9.5%) [O, 
Ho 67 ±3. 2 (4.8%) [O, 



,/i 2 prior] 
n h 2 prior] 



ory (|Crocce fe Scoccimarr j|2006l ; lMatsubarall2008T ) by 



fc, 



1 

6V 2 



dk Pn n (fc) 



-1/2 



(19) 



where Pn n (k) is again the linear power spectrum. ACDM 
predicts a value of fc, ~ 0.17/i Mpc . However, we will 
include fc* as a free fitting parameter. 

The scale dependance of the 6dFGS bias, B(s), is de- 
rived from the GiggleZ simulation (Poole et al., in prepara- 
tion); a dark matter simulation containing 2160 3 particles 
in a lh^ 1 Gpc box. We rank-order the halos of this simula- 
tion by V ma x and choose a contiguous set of 250 000 of them, 
selected to have the same clustering amplitude of 6dFGS as 
quantified by the separation scale ro, where £(ro) = 1. In 
the case of 6dFGS we found ro = 9.3/i _1 Mpc. Using the 
redshift space correlation function of these halos and of a 
randomly subsampled set of ~ 10 6 dark matter particles, 
we obtain 



B(s) = 1+ (s/0.474h -1 Mpc)" 



(20) 



which describes a 1.7% correction of the correlation function 
amplitude at separation scales of lO/t" 1 Mpc. To derive this 
function, the GiggleZ correlation function (snapshot 2 = 0) 
has been fitted down to Qh^ 1 Mpc, well below the smallest 
scales we are interested in. 



5 EXTRACTING THE BAO SIGNAL 

In this section we fit the model correlation function devel- 
oped in the previous section to our data. Such a fit can be 
used to derive the distance scale Dv(z e s) at the effective 
redshift of the survey. 



5.1 Fitting preparation 

The effective redshift of our sample is determined by 



(21) 



where Nb is the number of galaxies in a particular separa- 
tion bin and Wi and Wj are the weights for those galaxies 
from eq. El We choose z e a from bin 10 which has the limits 
lOO/i -1 Mpc and llO/i" 1 Mpc and which gave z cS = 0.106. 
Other bins show values very similar to this, with a standard 
deviation of ±0.001. The final result does not depend 
on a very precise determination of z e s, since we are not 
constraining a distance to the mean redshift, but a distance 
ratio (see equation 1241 later). In fact, if the fiducial model is 
correct, the result is completely independent of z e g. Only if 
there is a z-dependent deviation from the fiducial model do 
we need z e g to quantify this deviation at a specific redshift. 

Along the line-of-sight, the BAO signal directly con- 
strains the Hubble constant H(z) at redshift z. When mea- 
sured in a redshi ft shell, it constra ins the angular diameter 
distance Da(z) (|Matsubaral 1200-1 ) . In order to separately 
measure Da(z) and H(z) we require a BAO detection in 
the 2D correlation function, where it will appear as a ring 
at around 105k" 1 Mpc. Extremely large volumes are neces- 
sary for such a measurement. While there are studies that 
report a successful (but very low signal-to-noise) detection in 
the 2D correlation function using the SDSS-LRG data (e.g. 
iGaztanaga. Cabre fc Hul|2009l ; IChuang fc Wandl201il , but 
see also Kazin et al.l 201C ) , our sample does not allow this 
kind of analysis. Hence we restrict ourself to the ID corre- 
lation function, where we measure a combination of Da(z) 
and H(z). What we actually measure is a superposition of 
two angular measurements (R.A. and Dec.) and one line-of- 
sight measurement (redshift). To account for this mixture of 
measureme nts it is common to report the BAO distance con- 
strain ts as (jEisenstein et aL 1 120051 ; iPadmanabhan fc White! 
|200S| ) 



D v (z) 



(l + z) 2 D 2 A (z) 



H E(z) 



1/3 



(22) 



where Da is the angular distance, which in the case of fife = 
is given by Da(z) = Dc(z)/(1 + z). 

To derive model power spectra from CAMB we have to 
specify a complete cosmological model, which in the case of 
the simplest ACDM model (fife = 0, w = —1), is specified 
by six parameters: Lj c ,ujb,n a ,T, A a and h. These parame- 
ters are: the physical cold dark matter and baryon density, 
(uj c = fi c /i 2 , uib — fii,/i 2 ), the scalar spectral index, (n B ), the 
optical depth at recombination, (r), the scalar amplitude of 
the CMB temperature fluctuation, (As), and the Hubble 
constant in units of 100 km s~ 1 Mpc _1 (h). 

Our fit uses the parameter values from WMAP- 
7 (|Komatsu et all 120101 ): fi 6 fi 2 = 0.02227, r = 0.085 and 
Tis = 0.966 (maximum likelihood values). The scalar ampli- 
tude A s is set so that it results in erg = 0.8, which depends 
on fi m /i 2 . However erg is degenerated with the bias parame- 
ter 6 which is a free parameter in our fit. Furthermore, h is 
set to 0.7 in the fiducial model, but can vary freely in our 
fit through a scale distortion parameter a, which enters the 
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Figure 4. Likelihood contours of the distance Dy(z a g) against 
Q. m h 2 . The corresponding values of a are given on the right-hand 
axis. The contours show 1 and 2<r errors for both a full fit (blue 
solid contours) and a fit over 20 - 190/i _1 Mpc (black dashed 
contours) excluding the first data point. The black cross marks 
the best fitting values corresponding to the dashed black contours 
with (D v ,n m h 2 ) = (462,0.129), while the blue cross marks the 
best fitting values for the blue contours. The black solid curve 
corresponds to a constant Q m h 2 Dy (z e ff ) (Dy ~ /i -1 ), while the 
dashed line corresponds to a constant angular size of the sound 
horizon, as described in the text. 



model as 



£model(s) — £model ( as ) ■ 



(23) 



This parameter accounts for deviations from the fiducial cos- 
mological model, which we use to de rive distances from the 
measured redshift. It is defin ed as (|Eisenstein et all 120051 ; 
iPadmanabhan fc WhitelbOOST ) 



Dy{z c ff) 

Dy d (z cB )- 



(24) 



The parameter a enables us to fit the correlation function 
derived with the fiducial model, without the need to re- 
calculate the correlation function for every new cosmological 
parameter set. 

At low redshift we can approximate H(z) ~ Ho, which 
results in 



Ho 



(25) 



Compared to the correct equation [24] this approximation 
has an error of about 3% at redshift z = 0.1 for our fiducial 
model. Since this is a significant systematic bias, we do not 
use this approximation at any point in our analysis. 



5.2 Extracting Dy(z c a) and r a (z d )/ Dy(z c a) 

Using the model introduced above we performed fits to 18 
data points between 10/i -1 Mpc and 190/i _1 Mpc. We ex- 
cluded the data below 10/i -1 Mpc, since our model for non- 



linearities is not good enough to capture the effects on such 
scales. The upper limit is chosen to be well above the BAO 
scale, although the constraining contribution of the bins 
above 130/i -1 Mpc is very small. Our final model has 4 free 
parameters: Q, m h 2 , b, a and fc, . 

The best fit corresponds to a minimum \ 2 °f 15-7 with 
14 degrees of freedom (18 data-points and 4 free parame- 
ters) . The best fitting model is included in Figure [5] (black 
line). The parameter values are Qmh 2 = 0.138 ± 0.020, 
b = 1.81 ± 0.13 and a = 1.036 ± 0.062, where the er- 
rors are derived for each parameter by marginalising over 
all other parameters. For fc, we can give a lower limit of 
fc, = 0.19fe Mpc" 1 (with 95% confidence level). 

We can use eq. [24] to turn the measurement of a into 
a measurement of the distance to the effective redshift 
Dy(zeft) = otDy d (z c g) — 456 ± 27 Mpc, with a precision 
of 5.9%. Our fiducial model gives 4 d (z cff ) = 440. 5 Mpc 
where we have followed the distance definitions of I Wright! 
|2006f l throughout. For each fit we derive the parameter 
P = fl m {z) ' a45 /b, which we need to calculate the wide angle 
corrections for the correlation function. 

The maximum likelihood distribution of fc* seems to 
prefer smaller values than predicted by ACDM, although 
we are not able to constrain this parameter very well. This 
is connected to the high significance of the BAO peak in the 
6dFGS data (see Section[7]|. A smaller value of fc* damps the 
BAO peak and weakens the distance constraint. For compar- 
ison we also performed a fit fixing fc, to the ACDM predic- 
tion of fc* ~ 0.17h Mpc -1 . We found that the error on the 
distance Dy(z e a) increases from 5.9% to 8%. However since 
the data do not seem to support such a small value of fc, we 
prefer to marginalise over this parameter. 

The contours of Dy(z e ft) — fi m /i 2 are s hown in Figure [4] 
toget her with two degeneracy predictions (|Eisenstein et al] 
l2005h . The solid line is that of constant fi m /i 2 Dy(z e s), 
which gives the direction of degeneracy for a pure CDM 
model, where only the shape of the correlation function 
contributes to the fit, without a BAO peak. The dashed 
line corresponds to a constant r 3 (zd)/Dy(z e a), which is the 
degeneracy if only the position of the acoustic scale con- 
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tributes to the fit. The dashed contours exclude the first 
data point, fitting from 20 — 190/i _1 Mpc only, with the 
best fitting values a = 1.049 ± 0.071 (corresponding to 
D v (z cS ) = 462 ± 31 Mpc), Q m h 2 = 0.129 ± 0.025 and 
b = 1.72 ± 0.17. The contours of this fit are tilted towards 
the dashed line, which means that the fit is now driven by 
the BAO peak, while the general fit (solid contours) seems 
to have some contribution from the shape of the correlation 
function. Excluding the first data point increases the error 
on the distance constraint only slightly from 5.9% to 6.8%. 
The value of Q m h 2 tends to be smaller, but agrees within 
la with the former value. 

Going back to the complete fit from 10 — 190/t" 1 Mpc, 
we can include an external prior on il m h 2 from WMAP-7, 
which carries an error of only 4% (compared to the « 15% 
we obtain by fitting our data). Marginalising over Sl m h 2 now 
gives Dv{z c ff) — 459 ± 18 Mpc, which reduces the error 
from 5.9% to 3.9%. The uncertainty in Q m h 2 from WMAP- 
7 contributes only about 5% of the error in Dv (assuming no 
error in the WMAP-7 value of fl m h 2 results in Dv(z c g) = 
459 ± 17 Mpc). 



In Figure [5] we plot the ratio Dv(z)/D l 



'(*) 



function of redshift, where Dy~ z (z) = cz/H - At suffi- 
ciently low redshift the approximation H(z) ~ Ho is valid 
and the measurement is independent of any cosmological 
parameter except the Hubble constant. Th is figure also con- 
tains the results from lPercival et al.l (|2010l ). 

Rather than including the WMAP-7 prior on Q, m h 2 to 
break the degeneracy between Q m h 2 and the distance con- 
straint, we can fit the ratio r 3 (zd) /Dv(z e g), where r 3 (zd) is 
the sound horizon at the baryon drag epoch Zd- In princi- 
ple, this is rotating Figure 0] so that the dashed black line 
is parallel to the x-axis and hence breaks the degeneracy if 
the fit is driven by the BAO peak; it will be less efficient 
if the fit is driven by the shape of the correlation function. 
D uring the fit we calculat e r s (zd) using the fitting formula 
of lEisenstein fc Hul (119981 ). 

The best fit results in r s (z d )/ D v (z c(i ) = 0.336 ± 0.015, 
which has an error of 4.5%, smaller than the 5.9% found 
for Dv but larger than the error in Dy when adding the 
WMAP-7 prior on fl m h 2 . This is caused by the small dis- 
agreement in the Dv — il m h 2 degeneracy and the line of 
constant sound horizon in Figure [4] The \ 2 is 15.7, simi- 
lar to the previous fit with the same number of degrees of 
freedom. 



5.3 Extracting A{z c s) and i?(z e ff) 

We can also fit for the ratio of the distance between the 
effecti ve redshift, z e ff, and t he redshift of decoupling (z* = 
1091; lEisenstein et al.1120051 ); 



R(z c a) 



Dv(z e s) 



{i + z*)D A {z,y 



(26) 



with (1 + z,)Da(z*) being the CMB angular comoving dis- 
tance. Beside the fact that the Hubble constant Ho can- 
cels out in the determination of R, this ratio is also more 
robust against effects caused by possible extra relativistic 
species jEisenstein fc Whit jliool ) . We calculate D A {z,) for 
each fl m h 2 during the fit and then marginalise over Q m h 2 . 
The best fit results in R = 0.0324 ± 0.0015, with X 2 = 15.7 
and the same 14 degrees of freedom. 



Focusing on the path from z = to z e s = 0.106, our 
dataset can gi ve interesting constrain ts on Q m . We derive 
the parameter (|Eisenstein et al.1l2005l ) 



A(z eff ) = lOOLV^cff 



CZeff 



(27) 



which has no dependence on the Hubble constant since 
D v oc ft -1 . We obtain A(z cS ) = 0.526 ± 0.028 with 
X 2 /d.o.f. = 15.7/14. The value of A would be identical to 
vHm if measured at redshift z = 0. At redshift z e g = 0.106 
we obtain a deviation from this approximation of 6% for 
our fiducial model, which is small but systematic. We can 
express A, including the curvature term ilk and the dark 
energy equation of state parameter w, as 



A(z) = 



-, 2/3 



with 



and 



x {z) = Do{z)^ = 
c 




(28) 



dz 
E(z') 



e(z) = [n m (i + zf + n fc (i + z y 

+ H A (1 + 2 ) 3(1+ » ) )] 1/2 . 



(29) 



(30) 



Using this equation we now linearise our result for fi m in ilk 
and w and get 

Q m = 0.287 + 0.039(1 + w) + 0.039Q k ± 0.027. (31) 

For comparison, lEisenstein et al.l (|2005l ) found 

fl m = 0.273 + 0.123(1 + w) + 0.137n fc ± 0.025 (32) 

based on the SDSS LRG DR3 sample. This result shows the 
reduced sensitivity of the 6dFGS measurement to w and Q k ■ 



6 COSMOLOGICAL IMPLICATIONS 

In this section we compare our results to other studies and 
discuss the implications for constraints on cosmological pa- 
rameters. We first note that we do not see any excess cor- 
relation on large scales as found in the SDSS-LRG sample. 
Our correlation function is in agreement with a crossover to 
negative scales at lAQh~ l Mpc, as predicted from ACDM. 



6.1 Constraining the Hubble constant, Ho 

We now use the 6dFGS data to derive an estimate of 
the Hubble constant. We use the 6dFGS measurement of 
r 3 (z d )/D v (0.106) = 0.336 ± 0.015 and fit directly for the 
Hubble constant and Q m - We combine our measurement 
with a prior on Q m h 2 coming from the WMAP-7 Markov 
chain results (|Komatsu et al.l 20ld ) . Combining the cluster- 
ing measurement with Q, m h 2 from the CMB corresponds to 
the calibration of the standard ruler. 

We obtain values of Ho = 67±3.2 km s _1 Mpc _1 (which 
has an uncertainty of only 4.8%) and Q, m = 0.296 ± 0.028. 
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Table 2. uiCDM constraints from different datasets. Comparing the two columns shows the influence of the 6dFGS d ata point. The 
6dFG S data point reduces the error on w by 24% compared to WMAP-7+LRG which contains only the BAO data points of lPercival et all 
We assume flat priors of O.f 1 < f! m /i 2 < 0.16 and marginalise over fl m h 2 . The asterisks denote the free parameters in each fit. 



parameter 


WMAP-7+LRG 


WMAP-7±LRG±6dFGS 


Ho 


69.9±3.8(*) 


68.7 ± 2. 8(*) 




0.283 ± 0.033 


0.293 ± 0.027 


Q A 


0.717 ± 0.033 


0.707 ±0.027 


w 


-1.01 ± 0.17(*) 


-0.97±0.13(*) 



Table 3. Parameter constraints from WMAP7+BAO for (i) a flat ACDM model, (ii) an open ACDM (oACDM), (iii) a flat model 
with w = const. (uiCDM), and (iv) an open model with w = constant (ouiCDM). We assume flat priors of 0.11 < £lmh 2 < 0.16 and 
marginalise over Q m h 2 . The asterisks denote the free parameters in each fit. 



parameter 


ACDM 


oACDM 


ibCDM 


otuCDM 




69.2 ± 1.1(*) 


68.3 ± 1.7(*) 


68.7±2.8(*) 


70.4±4.3(*) 




0.288 ±0.011 


0.290 ±0.019 


0.293 ±0.027 


0.274 ± 0.035 


n k 


(0) 


-0.0036 ± 0.0060(*) 


(0) 


-0.013 ±0.010(*) 




0.712 ±0.011 


0.714 ±0.020 


0.707 ± 0.027 


0.726 ± 0.036 


w 


(-1) 


(-1) 


-0.97±0.13(*) 


-1.24 ± 0.39(*) 




0.15 0.2 0.25 0.3 0.35 0.4 



Figur e 6. The blu e con tours show the WMAP-7 Cl m h 2 
prior llKomatsu et al.l |2010|) . The black contour shows con- 
straints from 6dFGS derived by fitting to the measurement of 
rs(zd)/ Dy(z e s). The solid red contours show the combined con- 
straints resulting in Ho = 67 ± 3.2 km s — 1 Mpc — 1 and Q m = 
0.296 ±0.028. Combining the clustering measurement with f2 m /i 2 
from the CMB corresponds to the calibration of the standard 
ruler. 




Figure 7. The blue contours sh ows the WMAP-7 degeneracy in 
Ho and w llKomatsu et alj|2010h . highlighting the need for a sec- 
ond datasct to break the degeneracy. The black contours show 
constraints from BAO data inc o rpora ting the r 3 (zd) / Dy (z e g ) 
measurements of iPercival et al.l 1 120101 1 and 6dFGS. The solid 
red contours show the combined constraints resulting in w = 
—0.97 ± 0.13. Excluding the 6dFGS data point widens the con- 
straints to the dashed red line with w = —1.01 ± 0.17. 



Table [T] and Figure [5] summarise the results. The value of 
fl m agrees with the value we derived earlier (Section 15.311 . 

To combine our measurement with the latest CMB data 
we use the WMAP-7 distance priors, namely the acoustic 
scale 

/a = (1 + *.)22aC*>, (33) 
r s (z«) 



the shift parameter 

Ji= 100^^(1 +*.)£>a(*0 (34) 
c 

an d the redshift of d ecoupling z, (Tables 9 and 10 
in iKomatsu etHI |201Ch . This combined analysis reduces 
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the error further and yields Ho = 68.7 ± 1.5 km s 1 Mpc 1 
(2.2 %) and Q, m = 0. 29 ± .022 (7.6%). 

iPercival et al.l (|2010l ) determine a value of Hq = 
68.6 ± 2.2 km s ^Mpc" 1 usi n g SDS S-DRT, SDSS-LRG and 
2dFGRS, while iReid et all (|201Ch found H Q = 69.4 ± 
1.6 km s^Mpc" 1 using the SDSS-LRG sample and WMAP- 
5. In contrast to these results, 6dFGS is less affected by 
parameters like Q,k and w because of its lower redshift. In 
any case, our result of the Hubble constant agrees very 
well with earlier BAO analyses. Furthermore our result 
agrees with the la test CMB measureme nt of Ho = 70.3 ± 
2.5kms" 1 Mpc" 1 (iKomatsu et alj|20ich. 

The SHOES program (|Riess et al.ll201ll ) determined the 
Hubble constant using the distance ladder method. They 
used about 600 near-IR observations of Cepheids in eight 
galaxies to improve the calibration of 240 low redshift 
[z < 0.1) SN la, and calibrated the Cepheid distances us- 
ing the geometric distance to the maser galaxy NGC 4258. 
They found H = 73.8 ± 2.4 km s _1 Mpc _1 , a value con- 
sistent with the initial results of the Hubb le Ke y project 
iFreedman et all (Ho = 72±8 km s _1 Mpc _1 ; 1200 ll ) but 1. To- 
higher than our value (and 1.8cr higher when we combine 
our dataset with WMAP-T). While this could point toward 
unaccounted or under-estimated systematic errors in either 
one of the methods, the likelihood of such a deviation by 
chance is about 10% and hence is not enough to represent a 
significant discrepancy. Possible systematic errors affecting 
the BAO measurements are the modelling of non-linearities, 
bias and redshift-space distortions, although these system- 
atics are not expected to be significant at the large scales 
relevant to our analysis. 

To summarise the finding of this section we can state 
that our measurement of the Hubble constant is competitive 
with the latest result of the distance ladder method. The dif- 
ferent techniques employed to derive these results have very 
different potential systematic errors. Furthermore we found 
that BAO studies provide the most accurate measurement 
of Ho that exists, when combined with the CMB distance 
priors. 



6.2 Constraining dark energy 

One key problem driving current cosmology is the determi- 
nation of the dark energy equation of state parameter, w. 
When adding additional parameters like w to ACDM we find 
large degeneracies in the WMAP-T-only data. One example 
is shown in Figure [7) WMAP-T alone can not constrain Ho 
or w within sensible physical boundaries (e.g. w < — 1/3). As 
we are sensitive to Ho , we can break the degeneracy between 
w and Ho inherent in the CMB-only data. Our assumption 
of a fiducial cosmology with w = — 1 does not introduce a 
bias, since our data is not sensitive to this parameter and 
any deviation from this assumption is modelled within the 
shift parameter a. 

We again use the WMAP-T distance priors intro- 
duced in the last section. In addition to our value of 
r a ( z d )/D v (0.106) = .336 ± 0.015 we use the results 
of IPercival et ail (|2010l ). who found r a (z d )/D v (0.2) = 
0.1905 ± 0.0061 and r a (z d )/D v (0.35) = 0.109T ± 0.0036. 
To account for the correlation between the two latter data 
points we employ the covariance matrix reported in their 
paper. Our fit has 3 free parameters, Q m h 2 , Ho and w. 
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Figure 8. The number of log-normal realisations found with a 
certain \f Ax 2 , where the A.\ 2 is obtained by comparing a fit us- 
ing a ACDM correlation function model with a no-baryon model. 
The blue line indicates the 6dFGS result. 
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Figure 9. The different log-normal realisations used to calculate 
the covariance matrix (shown in grey). The red points indicate 
the mean values, while the blue points show actual 6dFGS data 
(the data point at 5h~ 1 Mpc is not included in the fit). The red 
data points arc shifted by 2h~ 1 Mpc to the right for clarity. 



The best fit gives w = -0.9T ± 0.13, Ho = 68. T ± 
2.8 kms^'Mpc^ 1 and ft m /i 2 = 0.1380 ± 0.0055, with a 
X 2 /d.o.f. = 1.3/3. Tableland Figure [7] summarise the re- 
sults. To illustrate the importance of the 6dFGS result to 
the overall fit we also show how the results change if 6dFGS 
is omitted. The 6dFGS data improve the constraint on w by 
24%. 

Finally we show the best fitting cosmological parame- 
ters for different cosmological models using WMAP-T and 
BAO results in Table [3] 



7 SIGNIFICANCE OF THE BAO DETECTION 

To test the signi ficance of our d e tectio n of the BAO sig- 
nature we follow lEisenstein et al.l (|2005l ) and perform a fit 
with a fixed fib ~ 0, which corresponds to a pure CDM 
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a 



Figure 10. This plot shows the distribution of the parameter a 
derived from the 200 log-normal realisations (black). The distri- 
bution is well fit by a Gaussian with a mean of fj, = 0.998 ± 0.004 
and a width of a = 0.057 ± 0.005. In blue we show the same dis- 
tribution selecting only the log-normal realisations with a strong 
BAO peak (> 2<r). The Gaussian distribution in this case gives a 
mean of 1.007 ± 0.007 and a = 0.041 ± 0.008. 



model without a BAO signature. The best fit has x 2 = 21.4 
with 14 degrees of freedom and is shown as the red dashed 
line in Figure [2] The parameter values of this fit depend on 
the parameter priors, which we set to 0.7 < a < 1.3 and 

0. 1 < Q. m h < 0.2. Values of a much further away from 
1 are problematic since eq. [24] is only valid for a close to 

1. Comparing the best pure CDM model with our previous 
fit, we estimate that the BAO signal is detected with a sig- 
nificance of 2.4cr (corresponding to A\ 2 = 5.6). As a more 
qualitative argument for the detection of the BAO signal we 
would like to refer again to Figure U where the direction of 
the degeneracy clearly indicates the sensitivity to the BAO 
peak. 

We can also use the log-normal realisations to deter- 
mine how likely it is to find a BAO detection in a survey 
like 6dFGS. To do this, we produced 200 log-normal mock 
catalogues and calculated the correlation function for each 
of them. We can now fit our correlation function model to 
these realisations. Furthermore, we fit a no-baryon model to 
the correlation function and calculate Ax , the distribution 
of which is shown in Figure [8] We find that 26% of all re- 
alisations have at feast a 2a BAO detection, and that 12% 
have a detection > 2.4a. The log-normal realisations show a 
mean significance of the BAO detection of 1.7±0.7ct, where 
the error describes the variance around the mean. 

Figure shows the 6dFGS data points together with 
all 200 log- normal realisations (grey). The red data points 
indicate the mean for each bin and the black line is the input 
model derived as explained in Section |3. 31 This comparison 
shows that the 6dFGS data contain a BAO peak slightly 
larger than expected in ACDM. 

The amplitude of the acoustic feature relative to 
the overall normalisation of the galaxy correlation func- 
tion is quite sensitive to the baryon fraction, fi = 
Qb/Qrn l|Matsubaral l2004D. A higher BAO peak could hence 
point towards a larger baryon fraction in the local universe. 
However since the correlation function model seems to agree 
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Figure 11. Redshift distribution of 6dFGS, WALLABY and two 
different versions of the proposed TAIPAN survey. See text for 
details. 



very well with the data (with a reduced \ 2 °f 1-12) and is 
within the range spanned by our log-normal realisations, we 
can not claim any discrepancy with ACDM. Therefore, the 
most likely explanation for the excess correlation in the BAO 
peak is sample variance. 

In Figure [10] we show the distribution of the parameter 
a obtained from the 200 log-normal realisations. The dis- 
tribution is well described by a Gaussian with x 2 /d.o.f. = 
14.2/20, where we employed Poisson errors for each bin. 
This confirms that a has Gaussian distributed errors in the 
approximation that the 6dFGS sample is well-described by 
log-normal realisations of an underlying ACDM power spec- 
trum. This result increases our confidence that the appli- 
cation of Gaussian errors for the cosmological parameter 
fits is correct. The mean of the Gaussian distribution is at 
0.998 ±0.004 in agreement with unity, which shows, that we 
are able to recover the input model. The width of the dis- 
tribution shows the mean expected error in a in a ACDM 
universe for a 6dFGS-like survey. We found a — 0.057±0.005 
which is in agreement with our error in a of 5.9%. Figure [TUI 
also contains the distribution of a, selecting only the log- 
normal realisations with a strong (> 2a) BAO peak (blue 
data). We included this selection to show, that a stronger 
BAO peak does not bias the estimate of a in any direc- 
tion. The Gaussian fit gives x 2 /d.o.f. = 5/11 with a mean of 
1.007 ± 0.007. The distribution of a shows a smaller spread 
with a = 0.041 ±0.008, about 2a below our error on a. This 
result shows, that a survey like 6dFGS is able to constrain 
q (and hence Dy and Ho) to the precision we report in this 
paper. 



8 FUTURE ALL SKY SURVEYS 

A major new wide-sky survey of the local Universe will be 
the Wide field ASKAP L-band Legacy All-sky Blind surveY 
(WALLABY0. This is a blind HI survey planned for the 
Australian SKA Pathfinder telescope (ASKAP), currently 

2 http: //www. atnf.csiro.au/research/WALLABY 
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Figure 12. Predictions for two versions of the proposed TAIPAN 
survey. Both predictions assume a 2n stcradian southern sky- 
coverage, excluding the Galactic plane (i.e. |6| > 10°). TAIPAN1 
contains 406 000 galaxies while TAIPAN2 contains 221000, (sec 
Figure ITTl l. The blue points are shifted by 2h~ 1 Mpc to the right 
for clarity. The black line is the input model, which is a ACDM 
model with a bias of 1.6, /3 = 0.3 and fc* = 0A7h Mpc -1 . For 
a large number of realisations, the difference between the input 
model and the mean (the data points) is only the convolution 
with the window function. 
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Figure 13. Prediction for the WALLABY survey. We have as- 
sumed a An stcradian survey with 602 000 galaxies, b = 0.7, 
P = 0.7 and fc* = 0.V7H Mpc" 1 . 



under construction at the Murchison Radio-astronomy Ob- 
servatory (MRO) in Western Australia. 

The survey will cover at least 75% of the sky with the 
potential to cover 4tt of sky if the Westerbork Radio Tele- 
scope delivers complementary northern coverage. Compared 
to 6dFGS, WALLABY will more than double the sky cov- 
erage including the Galactic plane. WALLABY will contain 
~ 500 000 to 600 000 galaxies with a mean redshift of around 
0.04, giving it around 4 times greater galaxy density com- 
pared to 6dFGS. In the calculations that follow, we assume 
for WALLABY a 4-7T survey without any exclusion around 
the Galactic plane. The effective volume in this case turns 
out to be 0.12/i" 3 Gpc 3 . 



The TAIPAN surve£| proposed for the UK Schmidt 
Telescope at Siding Spring Observatory, will cover a com- 
parable area of sky, and will extend 6dFGS in both depth 
and redshift (z ~ 0.08). 

The redshift distribution of both surveys is shown in 
Figure [TT1 alongside 6dFGS. Since the TAIPAN survey is 
still in the early planning stage we consider two realisa- 
tions: TAIPAN1 (406 000 galaxies to a faint magnitude limit 
of r = 17) and the shallower TAIPAN2 (221000 galax- 
ies to r = 16.5). We have adopted the same survey win- 
dow as was used for 6dFGS, meaning that it covers the 
whole southern sky excluding a 10° strip around the Galac- 
tic plane. The effective volumes of TAIPAN 1 and TAIPAN 2 
are 0.23/i~ 3 Gpc 3 and 0.13/i~ 3 Gpc 3 , respectively. 

To predict the ability of these surveys to measure the 
large scale correlation function we produced 100 log-normal 
realisations for TAIPAN 1 and WALLABY and 200 log- 
normal realisations for TAIPAN2. Figures [12] and [13] show 
the results in each case. The data points are the mean of 
the different realisations, and the error bars are the diago- 
nal of the covariance matrix. The black line represents the 
input model which is a ACDM prediction convolved with 
a Gaussian damping term using k, = 0.17/1 Mpc- 1 (see 
eq. I17|) . We used a bias parameter of 1.6 for TAIPAN and 
following our fiducial model we get f3 = 0.3, resulting in 
A = b 2 (l + 2/3/3 + /3 2 /5) = 3.1. For WALLABY we used a 
bias o f 0.7 (based on the r esults found in the HIPASS sur- 
vey; iBasilakos et alj|2007h . This results in f) = 0.7 and 
A = 0.76. To calculate the correlation function we used 
Po = 40 000/i 3 Mpc 3 for TAIPAN and P = 5 000/i 3 Mpc 3 
for WALLABY. 

The error bar for TAIPAN 1 is smaller by roughly a fac- 
tor of 1.7 relative to 6dFGS, which is consistent with scal- 
ing by y/Vcii and is comparable to the SDSS-LRG sample. 
We calculate the significance of the BAO detection for each 
log-normal realisation by performing fits to the correlation 
function using ACDM parameters and — 0, in exactly 
the same manner as the 6dFGS analysis described earlier. 
We find a 3.5 ± 0.8<r significance for the BAO detection for 
TAIPAN 1, 2.1±0.7ff for TAIPAN2 and 2.1±0.7cr for WAL- 
LABY, where the error again describes the variance around 
the mean. 

We then fit a correlation function model to the mean 
values of the log-normal realisations for each survey, using 
the covariance matrix derived from these log-normal real- 
isations. We evaluated the correlation function of WAL- 
LABY, TAIPAN 2 and TAIPAN 1 at the effective redshifts 
of 0.1, 0.12 and 0.14, respectively. With these in hand, we 
are able to derive distance constraints to respective preci- 
sions of 7%, 6% and 3%. The predicted value for WAL- 
LABY is not significantly better than that from 6dFGS. 
This is due to the significance of the 6dFGS BAO peak in 
the data, allowing us to place tight constraints on the dis- 
tance. As an alternative figure-of-merit, we derive the con- 
straints on the Hubble constant. All surveys recover the in- 
put parameter of Ho = 70 km s Mpc , with absolute un- 
certainties of 3.7, 3 and 2.2 kms'^pc" 1 for WALLABY, 
TAIPAN2 and TAIPAN 1, respectively. Hence, TAIPAN 1 is 



3 TAIPAN: Transforming Astronomical Imaging surveys through 
Polychromatic Analysis of Nebulae 
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able to constrain the Hubble constant to 3% precision. These 
constraints might improve when combined with Planck con- 
straints on f2i,/i 2 and Q, m h 2 which will be available when 
these surveys come along. 

Since there is significant overlap between the survey 
volume of 6dFGS, TAIPAN and WALLABY, it might be 
interesting to test whether the BAO analysis of the local 
universe can mak e use of a multiple tracer analysis, as sug- 
gested recently bv lArnalte-Mur et al.l (|201lh . These authors 
claim that by employing two different tracers of the matter 
density field - one with high bias to trace the central over- 
densities, and one with low bias to trace the small density 
fluctuations - one can i mprove the detection and measure- 
ment of the BAO signal. lArnalte-Mur et al.l (|201lh test this 
approach using the SDSS-LRG sample (with a very large 
bias) and the SDSS-main sample (with a low bias). Although 
the volume is limited by the amount of sample overlap, they 
detect the BAO peak at 4.1a. Likewise, we expect that the 
contrasting high bias of 6dFGS and TAIPAN, when used in 
conjunction with the low bias of WALLABY, would furnish 
a combined sample that would be ideal for such an analysis. 

Neither TAIPAN nor WALLABY are designed as BAO 
surveys, with their primary goals relating to galaxy forma- 
tion and the local universe. However, we have found that 
TAIPAN 1 would be able to improve the measurement of the 
local Hubble constant by about 30% compared to 6dFGS go- 
ing to only slightly higher redshift. WALLABY could make 
some interesting contributions in the form of a multiple 
tracer analysis. 



9 CONCLUSION 

We have calculated the large-scale correlation function of the 
6dF Galaxy Survey and detected a BAO peak with a sig- 
nificance of 2.4<r. Although 6dFGS was never designed as a 
BAO survey, the peak is detectable because the survey con- 
tains a large number of very bright, highly biased galaxies, 
within a sufficiently large effective volume of 0.08/i~ 3 Gpc 3 . 
We draw the following conclusions from our work: 

• The 6dFGS BAO detection confirms the finding by 
SDSS and 2dFGRS of a peak in the correlation function at 
around 105/i _1 Mpc, consistent with ACDM. This is impor- 
tant because 6dFGS is an independent sample, with a dif- 
ferent target selection, redshift distribution, and bias com- 
pared to previous studies. The 6dFGS BAO measurement is 
the lowest redshift BAO measurement ever made. 

• We do not see any excess correlation at large scales 
as seen in the SDSS-LRG sample. Our correlation func- 
tion is consistent with a crossover to negative values at 
140/i _1 Mpc, as expected from ACDM models. 

• We derive the distance to the effective redshift as 
D v (z c a) = 456 ± 27 Mpc (5.9% precision). Alternatively, 
we can derive r a (zd) /Dy{z e s) = 0.336 ± 0.015 (4.5% preci- 
sion). All parameter constraints are summarised in Table [T] 

• Using a prior on Q m h 2 from WMAP-7, we find Q, m = 
0.296 ± 0.028. Independent of WMAP-7, and taking into 
account curvature and the dark energy equation of state, we 
derive Q m = 0.287 + 0.039(1 + w) + 0.039fi fc ± 0.027. This 
agrees very well with the first value, and shows the very 
small dependence on cosmology for parameter derivations 
from 6dFGS given its low redshift. 



• We are able to measure the Hubble constant, Ho = 
67 ± 3.2 km s _1 Mpc _1 , to 4.8% precision, using only the 
standard ruler calibration by the CMB (in form of 0. m h 2 and 
Qbh 2 ). Compared to previous BAO measurements, 6dFGS is 
almost completely independent of cosmological parameters 
(e.g. Qk and w), similar to Cepheid and low- z supernovae 
methods. However, in contrast to these methods, the BAO 
derivation of the Hubble constant depends on very basic 
early universe physics and avoids possible systematic errors 
coming from the build up of a distance ladder. 

• By combining the 6dFGS BAO measure ment with those 
of WMAP-7 and previous redshift samples Percival et al.l 
(from SDSS-DR7, SDDS-LRG and 2dFGRS; 120101 ), we can 
further improve the constraints on the dark energy equation 
of state, w, by breaking the Ho — w degeneracy in the CMB 
data. Doing this, we find w = —0.97 ± 0.13, which is an 
improvement of 24% compared to previous combinations of 
BAO and WMAP-7 data. 

• We have made detailed predictions for two next- 
generation low redshift surveys, WALLABY and TAIPAN. 
Using our 6dFGS result, we predict that both surveys will 
detect the BAO signal, and that WALLABY may be the first 
radio galaxy survey to do so. Furthermore, we predict that 
TAIPAN has the potential to constrain the Hubble constant 
to a precision of 3% improving the 6dFGS measurement by 
30%. 
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APPENDIX A: GENERATING LOG-NORMAL 
MOCK CATALOGUES 

Here we explain in detail the different steps used to de- 
rive a log-normal mock catalogue, as a useful guide for re- 
searchers in the field. We start with an input power spec- 
trum, (which is determined as explained in Section 13. 3p in 
units of ft~ 3 Mpc 3 . We set up a 3D grid with the dimensions 
L x x L y x L z = 1000 x 1000 x 1000ft -1 Mpc with 200 3 sub- 
cells. We then distribute the quantity P(k) /V over this grid, 
where V is the volume of the grid and k = yfk% + ky + fc| 
with k x = n x 2n/L x and n x being an integer value specifying 
the x coordinates of the grid cells. 

Performing a complex-to-real Fourier transform (FT) of 
this grid will produce a 3D correlation function. Since the 
power spectrum has the property P(—k) — P(k)* the result 
will be real. 

The next step is to replace the correlation function £(r) 
at each point in the 3D grid by ln[l + £(r)], where In is 
the natural logarithm. This step prepares the input model 
for the inverse step, which we later use to produce the log- 
normal density field. 

Using a real-to-complex FT we can revert to fc-space 
where we now have a modified power spectrum, P\ n (k). 
At this point we divide by the number of sub-cells N c . 
The precise normalisation depends on the definition of 
the di screte Fourier transfo rm. We use the FFTW li- 
brary (|Frigo fc Johnsonll2005l ). where the discrete FT is de- 
fined as 

N c -1 

Yi = ]T Xj exp [±2irij yf^l/N c ] . (Al) 

3=0 

The modified power spectrum Pi n (A?) is not guarantied to be 
neither positive defined nor a real functi on, which contra- 
dicts the definition of a power spectrum. IWeinberg fc Cold 
(|1992| ) suggested to construct a well defined power spectrum 
from Pin(fc) by 

P{ n (k) = max [o, Re[P ln (fc)]] . (A2) 

We now generate a real and an imaginary Fourier amplitude 
S(k) for each point on the grid by randomly sampling from 

a Gaussian distribution with r.m.s. \J P/ n (fc)/2. However, to 
ensure that the final over-density field is real, we have to 
manipulate the grid, so that all sub-cells follow the condition 
5{-%) = 5(h)*. 

Performing another FT results in an over-density field 
S(x) from which we calculate the variance a%. The mean 
of 8(x) should be zero. The log- normal density field is then 
given by 

^ (x) = exp [8(x)- a%/2] , (A3) 




20 40 60 80 100 120 140 160 180 200 
s [h" 1 Mpc] 

Figure Bl. Correlation function error for different values of P$. 
The weighting with Po = 40 000/i 3 Mpc -3 reduces the error at 
the BAO scale by almost a factor of four compared to the case 
without weighting. The red dashed line indicates the jack-knife 
error. 



which is now a quantity defined on [0, oo[ only, while S(x) is 
defined on ] — oo, oo[. 

Since we want to calculate a mock catalogue for a par- 
ticular survey we have to incorporate the survey selection 
function. If W(x) is the selection function with the normali- 
sation ^ W(x) — 1, we calculate the mean number of galax- 
ies in each grid cell as 

n g {g) = NW{$) (A4) 

where N is the total number of galaxies in our sample. The 
galaxy catalogue itself is than generated by Poisson sampling 

Tlg(x). 

The galaxy position is not defined within the sub-cell, 
and we place the galaxy in a random position within the box. 
This means that the correlation function calculated from 
such a distribution is smooth at scales smaller than the sub- 
cell. It is therefore important to make sure that the grid 
cells are smaller than the size of the bins in the correlation 
function calculation. In the 6dFGS calculations presented in 
this paper the grid cells have a size of 5/i _1 Mpc, while the 
correlation function bins are Wh- 1 Mpc in size. 



APPENDIX B: COMPARISON OF 
LOG-NORMAL AND JACK-KNIFE ERROR 
ESTIMATES 

We have also estimated jack-knife errors for the correlation 
function, by way of comparison. We divided the survey into 
18 regions and calculated the correlation function by ex- 
cluding one region at a time. We found that the size of the 
error-bars around the BAO peak varies by around 20% in 
some bins, when we increase the number of jack-knife regions 
from 18 to 32. Furthermore the covariance matrix derived 
from jack-knife resampling is very noisy and hard to invert. 

We show the jack-knife errors in Figure IB1I The jack- 
knife error shows more noise and is larger in most bins com- 
pared to the log-normal error. The error shown in Figure iBll 
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APPENDIX C: WIDE-ANGLE FORMALISM 
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Figure B2. Correlation matrix of the jack-knife errors (upper 
left triangle) and log-normal errors (lower right triangle). 



is only the diagonal term of the covariance matrix and does 
not include any correlation between bins. 

The full error matrix is shown in Figure IB2I where we 
plot the correlation matrix of the jack-knife error estimate 
compared to the log-normal error. The jack-knife correla- 
tion matrix looks much more noisy and seems to have less 
correlation in neighbouring bins. 

The number of jack-knife regions can not be chosen ar- 
bitrarily. Each jack-knife region must be at least as big as 
the maximum scale under investigation. Since we want to 
test scales up to almost 200/i~ 1 Mpc our jack-knife regions 
must be very large. On the other hand we need at least as 
many jack-knife regions as we have bins in our correlation 
function, otherwise the covariance matrix is singular. These 
requirements can contradict each other, especially if large 
scales are analysed. Furthermore the small number of jack- 
knife regions is the main source o f noise (for a more de tailed 
study of jack-knife errors see e.g. iNorberg et al"1l2008h . 

Given these limitations in the jack-knife error approach, 
correlation function studies on large scales usually employ 
simulations or log-normal realisations to derive the covari- 
ance matrix. We decided to use the log-normal error in our 
analysis. We showed that the jack-knife errors tend to be 
larger than the log-normal error at larger scales and carry 
less correlation. These differences might be connected to the 
much higher noise level in the jack-knife errors, which is 
clearly visible in all our data. It could be, however, that 
our jack-knife regions are too small to deliver reliable errors 
on large scales. We use the minimum number of jack-knife 
regions to make the covariance matrix non-singular (the cor- 
relation function is measured in 18 bins). The mean distance 
of the jack-knife regions to each other is about 200/i _1 Mpc 
at the mean redshift of the survey, but smaller at low red- 
shift. 



The general redshift space correlation function (ignoring the 
plane parallel approximation) depends on </>, 8 and s. Here, 
s is the separation between the galaxy pair, 8 is the half 
opening angle, a nd 4> is the angle of s to the line of sight 
(see Figure 1 in iRaccanelli et alj|20ich . For the following 
calculations it must be considered that in this parametrisa- 
tion, (j> an d 6 are not independent. 

The total correlation functio n model, including Q(8 2 ) 
correction terms, is then given bv lPapai fc Szapudil (|2008l ). 



6 



, s) = aoo + 2ao2 cos(20) + 022 cos(20) + 622 sin (20) 
+ ^ — 4a 2 cos(20) — 4a22 — 4&22 — 4aio cot 2 (0) 



+ 4au cot (0) - 4ai2 cot (0) cos(20) + 4b 



11 



b 12 cos 2 (</>)] 8 2 + 0{8 i ) 



(CI) 



This equation reduces to the plane parallel approximation if 
8 = 0. The factors a xy and b xy in this equation are given by 

• 2 - 



aoo 



a,02 = 



, 2/3 2/3^ 
1 + — + — 
3 15 



£o(r) 



p 2/P 
3 21 



t-i, \ , 3/3 2 2/ \ 
6(r)+ 140 



P_ W_ 
2 14 



«22 



&22 = 



aio 



t 

15 



e (r) 



p 2 



21 6(r) + W^ (r) 



t 
21 



t-2ss 4/3 2/ N 

6(>) - - 35-61 M 



2/3 + 



4/3 2 



(C2) 



5r 



an 



«21 = 



bit = 



&12 = 



|Uo°(r)-26V) 
fJ[3C3 1 M-2C 1 1 (r) 



where P = fl m (z) ' 545 /b, with b being the linear bias. The 
correlation function moments are given by 



= JL / dk k n IU(k)ji(rk) 



(C3) 



with ji(x) being the spherical Bessel function of order I. 

The final spherically averaged correlation function is 
given by 

C( S ) = r r' 2 a<P, 0, s)N(4>, e, s) Md4, (C4) 



where the function N(<j>, 8, s) is obtained from the data. 
N(<j>, 9, s) counts the number of galaxy pairs at different cj>, 8 
and s and includes the areal weighting sin(0) which usually 
has to be included in an integral over <f>. It is normalised 
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s [lr 1 Mpc] 

Figure CI. The half opening angle 8 as a function of separation 
s of the 6dFGS weighted catalogue. The plane parallel approxi- 
mation assumes 9 = 0. The mean half opening angle at the BAO 
scale is < 10°. The colour bar gives the number of pairs in each 
bin. 




l 
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Figure C2. The black line represents the plain correlation func- 
tion without redshift space distortions (RSD), £(r), obtained by 
a Hankel transform of our fiducial ACDM power spectrum. The 
blue line includes the linear model for redshift space distortions 
(linear Kaiser factor) using /3 = 0.27. The red line uses the same 
value of /3 but includes all correction terms outlined in eq. IC1I 
using the N(<j>, 9, s) distribution of the weighted 6dFGS sample 
employed in this analysis. 



such that 

/ / N{4>, 0, s) d6d<f> = 1. (C5) 
Jo Jo 

If the angle 6 is of order 1 rad, higher order terms become 
dominant and eg. I CI I is no longer sufficient. Our weighted 
sample has only small values of 8, but growing with s (see 
figure [CT|) . In our case the correction terms contribute only 
mildly at the BAO scale (red line in figure [C2| . However 
these corrections behave like a scale dependent bias and 
hence can introduce systematic errors if not modelled cor- 
rectly. 
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